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Abstract 

The adsorption of a single multi-block AS-copolymer on a solid planar substrate is investi- 
gated by means of computer simulations and scaling analysis. It is shown that the problem can 
be mapped onto an effective homopolymer adsorption problem. In particular we discuss how 
the critical adsorption energy and the fraction of adsorbed monomers depend on the block 
length M of sticking monomers A, and on the total length N of the polymer chains. Also 
the adsorption of the random copolymers is considered and found to be well described within 
the framework of the annealed approximation. For a better test of our theoretical prediction, 
two different Monte Carlo (MC) simulation methods were employed: a) off-lattice dynamic 
bead-spring model, based on the standard Metropolis algorithm (MA), and b) coarse-grained 
lattice model using the Pruned-enriched Rosenbluth method (PERM) which enables tests for 
very long chains. The findings of both methods are fully consistent and in good agreement 
with theoretical predictions. 

1 Introduction 

Adsorption of polymers on surfaces plays a key role in numerous technological applications and is 



^ also relevant to many biological processes. During the last three decades it has been constantly a 
5^ , focus of research interest. The theoretical studies of the behavior of polymers interacting with solid 
substrate have been based predominantly on both scaling analysis^"^ as well as on the self-consistent 
field (SCF) approach.'' The close relationship between theory and computer experiments in this 
field^'^ has proved especially fruitful. Most investigations focus as a rule on the determination of the 
critical adsorption point (CAP) location and on the scaling behavior of a variety of quantities below, 
above and at the CAP. Thus an eminent relation between polymer statistics and the corresponding 
correlation functions^ in the n-vector model of magnets with a free surface in the limit n — ^ has 
lead to a number of important results. Special interest has been payed to the determination of 
the so called crossover exponent (p which is known to govern the fraction of adsorbed monomers at 
the CAP. Recently the scaling relationship for a single chain adsorption has been tested by Monte 



Carlo (MC) simulation on a cubic lattice^'^ as well as by an off-lattice model^'^'^ and the adsorption 
transition of a polymer could be viewed nowadays as comparatively well understood. 

While the investigations mentioned above have been devoted exclusively to homopolymers, the 
adsorption of copolymers (e.g. multi-blocks or random copolymers) is still much less understood. 
Thus, for instance, the CAP dependence on block size M at fixed concentration of the sticking ti- 
mers is still unknown as are the scaling properties of regular multi-block copolymers in the vicinity of 
the CAP. Prom the theoretical perspective, the case of diblock copolymers has been studied mainly 
within the SCF-approach.^'^^ The case of random copolymers adsorption has gained comparatively 
more attention by researcher so far. It has been investigated by Whittington et al.^^'^^ using both 
the annealed and quenched models of randomness. In the latter case the authors implemented 
the Morita approximation (which is reduced to an optimization problem with a set of constraints 
involving the moments of the quenched random probability distribution) . The influence of sequence 
correlations on the adsorption of random copolymers has been studied by means of the variational 
and replica method approach.^'* Sumithra and Baumgaertner^^ examined the question of how the 
critical behavior of random copolymers differs from that of homopolymers. Thus, among a number 
of important conclusions, the results of Monte Carlo simulations demonstrated that the crossover 
exponent (see below) is independent of the fraction of attractive monomers /. 

In the present paper we use scaling analysis as well as two MC-simulation methods to study the 
critical behavior of multi-block and random copolymers. It turns out that the critical behaviour of 
these two types of copolymers could be reduced to the behavior of an effective homopolymer chain 
with " renormalized" segments. Por the multi- block copolymer this allows e.g. to explain how the 
critical attraction energy depends on the block length M and to derive an adsorption phase diagram 
in terms of CAP against M. In the case of random copolymers the sequence of sticky and neutral 
(as regards the solid substrate) monomers within a particular chain is fixed which exemplifies a 
system with quenched randomness. Nevertheless, close to criticality the chain is still rather mobile, 
so that the sequence dependence is effectively averaged over the time of the experiment and the 
problem can be reduced to the case of annealed randomness. We show that our MC-findings close 
to criticality could be perfectly treated within the annealed randomness model. 



2 Scaling properties of homopolymer adsorption 
2.1 Order parameter 

Before discussing copolymers adsorption we briefly sketch the scaling theory of homopolymer ad- 
sorption. It is well known that a single polymer chain undergoes a transition from a non-bound 
into an adsorbed state when the adsorption energy e per monomer increases beyond a critical value 
Cc- Here and in what follows e is measured in units of the thermal energy ksT (with Ub being the 
Boltzmann constant, ane T - the temperature of the system). The adsorption transition can be 
interpreted as a second-order phase transition at the critical point (CAP) of adsorption e = ec in the 
thermodynamical limit, i.e. iV — > oo. Close to the CAP the number of surface contacts Ng scales as 
Ns{e — €c) ~ N'^. The numerical value of (f) is somewhat controversial and hes in a range between 
(f) — 0.59 (ref.^) and = 0.484 (ref.^), we adopt however the value 4> = 0.50 ± 0.02 which has been 
suggested as the most satisfactory^'' by comparison with comprehensive simulation results. 

Consider a chain tethered to the surface at the one end. The fraction of monomers on the 
surface / = Ng/N may be viewed as an order parameter measuring the degree of adsorption. In the 
thermodynamic limit N — > oo, the fraction / goes to zero (~ 0{1/N)) for e << Sc, then near Sc, 
f ~ N'l'~^, and for e » ec (in the strong couphng hmit) / it is independent of N. Let us measure 
the distance from the CAP by the dimensionless quantity k — (e — ec)/ ec and also introduce the 
scaling variable rj = kN'^. The corresponding scaling ansatz is then 



fiv) = Nf-' G iv) . 



(1) 



with the scaling function 




for 77 — > 
for 77 » 1 



(2) 



The resulting scaling behavior of / follows as. 




for K «0 
for K, ^ 
for K » 1 



(3) 



2.2 Gyration radius 

The gyration radius in direction perpendicular to the surface, Rg±{r]), has the form 

R,±(V) = aN'^g,^ (rj) (4) 

One may determine the form of the scahng function Gg±{'r)) from the following consideration. At 
K < one has Rg± ~ aN'^, so that Qg± — const. In the opposite limit r] the A'"-dependence 
drops out and Qg±{ri) ~ ri~'^^^. Thus 

/const , for 77 < 
= |,-^/. , for ,»0 



As a result 



s±v /y ^_,/^ , for > ^ ^ 



The gyration radius in direction parallel to the surface has similar scaling representation: 

Rg\\{ri)^aN^gg\\{7^) (7) 

Again at k < the gyration radius Rg\\ ~ aN'^ and Qg\\ = const. At 77 ^ the chain behaves as 
a two-dimensional self-avoiding walk (SAW), i.e. Rg\\ ~ aN'"^, where 1^2 = 3/4 denotes the Flory 
exponent in two dimensions. In result, the scaling function behaves as 



^.11 (^) 

Thus 



const , at 77 < 

^(u2-u)/<t> ^ at 77>0 



„ , , 1 aN" , at ?7 < 

911 v/y {i^2-u)/4>j^u2 , at 77>0 



2.2.1 Blob picture 



In the limit kN'^ ^ 1 the adsorbed chain can be visualized as a string of adsorption blobs which 
forms a pancake-like quasi-two-dimensional layer on the surface. The blobs are defined to contain 
as many monomers g as necessary to be on the verge of being adsorbed and therefore carry an 
adsorption energy of the order of ksT each. The thickness of the pancake Ro\ corresponds to 



the size of the blob and the chain conformation within a blob stays unperturbed (i.e. it is simply 
a SAW), thus g ~ iRg±/a)^^'^ = k,^^^'^ where we have used eq [H The gyration radius can be 
represented thus as 

/ ]\r\ 

Rg\\ =Rg±{ — ) OC K^'^^-^^/^N"^ (10) 



9 

and one goes back to eq M which proves the consistency of the adsorption blob picture. Generally 
speaking, the number of blobs, N/ g ~ k^^'^N, is essential for the main scaling argument in the above- 
mentioned scaling functions. For example we could recast the order parameter scaling behavior eq[T] 
as 

/ = N^-'H (-] (11) 



9 



where H{x) denotes a new scaling function : 



I X , for X 1 

2.2.2 Ratio of gyration radius components 

The study of the ratio, r{r]) = Rg±/Rg\\, of gyration radius components is a convenient way to find 
the value of ec (see^'^°). In fact, from the previous scaling equations 

Hence at the critical point, i.e. at 77 — 0, the ratio r(0) = const is independent of A^. Thus by 
plotting r vs. e for different N all such curves should intersect at a single point which gives ec- 

Another way to fix ec is the following. Exactly at the critical point / ~ A^*^"^, so that by plotting 
fN^~'^ vs. A^ at different values of e one can determine the value e ~ ec under which fN^~'f' becomes 
independent of A^. 

2.3 Free energy of adsorption 

The adsorption on a surface at k > is due to a free energy gain which is proportional to the 
number of blobs, i.e., 



bulk 



AT 



OC K^/'^ . (14) 



The expression for the specific heat per monomer follows immediately from eq [H] as 



Cv 



{F — -Fbuik) 



—a 



(15) 



where a = 2 — (p^^. Note that a factor of fc^T is absorbed in the free energy throughout the paper. 
If = 0.5 then a = and the specific heat undergoes a jump at the CAP (cf. Section [6.1.2p . 

For a chain (of the length A^) on the verge of adsorption, the foregoing free energy gain, F — Fbuik, 
should be of the order of unity. In view of eq [T3] this gives an estimate for the critical energy of 
adsorption - CAP, 



where we have explicitly marked the CAP, ec{N) and ec(oo), for finite and infinitely long chains 
respectively. 

3 Multi-block copolymer adsorption 

Consider now the adsorption of a regular multi-block copolymer which is built up from monomers 
A which attract (stick) to the substrate and monomers B which are neutral to the substrate. In 
order to treat the adsorption of a regular multi-block AB - copolymer we reduce the problem to 
that of a homopolymer which has been considered above. The idea is that a regular multi-block 
copolymer can be considered as a "homopolymer" where a single AS-diblock plays the role of an 
effective monomer. ^'^ For such a mapping we first estimate the effective energy of adsorption per 
diblock. 

3.1 Effective energy of adsorption per diblock 

Each individual diblock is made up of an attractive A-block of length M and a neutral i?-block of 
the same length M. Upon adsorption the attractive A-block forms a string of blobs whereas the 
i?-part forms a non-adsorbed tail (or loop) - (see Figured]). The free energy gain of the attractive 
block may be written according to eq [TH as 




(16) 



attr 



(17) 



Figure 1: Schematic representation of an individual adsorbed Ai?-diblock. The A - part forms a 
string of quasi-two dimensional blobs and the B-pait is neutral regarding the substrate and its 
contribution to the free energy is of pure entropical nature. 

where we measure the energy in units of fc^T and n = (e — e^)/ measures the normalized distance 
from the CAP ej? of a homopolymer. The neutral 5-part which is most frequently a loop connecting 
adjacent ^-blocks, but could also be a tail with the one end free, contributes only to the entropy 
loss 

i"rep= (7-711) In M (18) 

where the universal exponents 7 and 711 are well known^'' (e.g. in 3D - space 7 = 1.159, 711 = 
—0.390). In case that also the tails are involved, one should also use the exponent 71 = 0.679 
albeit this does not change qualitatively the expression eq [181 They enter the partition function 
expressions for a free chain, a chain with both ends fixed at a two points, and for a chain, tethered 
by the one end.^^ In result the effective adsorption energy of a diblock is 

E{M) = K^/'*'M - (7 - 7ii) In M (19) 

3.2 Order parameter 

Now we consider a 'homopolymer' which is build up from effective units (diblocks), with the attrac- 
tive energy given by eqUHl Let us denote the total number of such effective units by A/" = N/2M. 
The fraction of effective units on the surface obeys then the same scaling law as given by eqlH i.e., 

^=U't>-^G{dJ\f'^) (20) 

where now ^ = [E — E^)/E^ with the critical adsoption energy E^ of the renormalized homopoly- 
mer. Generallv, one would expect E!^ to be of the order of e'i albeit for different models both critical 



energies would probably differ from each other. Eq [20] is accurate if one require that (i) k ^ 1 but 
M > 1 such that InM > 1 and k^/'^M > 1, and (ii) A/" > 1. The effective attraction of a 
segment of the renormailzed chain now depends on M according to eq [13 

Within each effective unit only Mg A-monomers will be adsorbed at criticality whereby this 
monomer number scales as 

M, = M'^G (kM^) (21) 

with K = (e - e^)/ej?. 

The total number of adsorbed monomer is given by 

Ns = MsMs = MsM^G [kW^) (22) 

It follows that the fraction 



- iM-'G(.M*)(|^)"'G(A(^)*) , ,23) 

where we have used the scaling law, eq [201 for the effective units. Hence, the final expression for 
the order parameter can be written as follows: 



/^liV*-.G(.M*)G(A(0) 



(24) 



Thus we have expressed the order parameter / of a multi-block copolymer in terms of the chain 
length A^, the block length M, the monomer attraction energy e as well as the model-dependent 
homopolymer critical attraction energy e^. Let us consider now some limiting cases. 

3.2.1 Close to criticality A = 

At the CAP of the multiblock chain one has A = 0, thus one can estimate the deviation k^, of the 
corresponding critical energy of adsorption, e^^, from that of a homopolymer, namely 

M _ -4 _ / (7-7ii)lnM + En 

where we have used eq [T^ and set (p = 0.5. Under this condition the second G -function in eq[21|is 
a constant, i.e., G(0) = const. On the other hand, with respect to a single effective unit the chain 



stays far from the criticality because of /tf = v/(7-7ii)lnM + E^^ > 1. In this case the first 
G - function in eq [24] behaves as G{k^^ \fM) ~ k^^ \/M where now is fixed by eq[25l In result, 
eq m] becomes 

;^{kl^Jl,}^^E±^f ,26) 



N 

3.2.2 State of the strong adsorption 

In this regime k^M > 1 and A^/N/M > 1 so that / ^ {l/y/N)G{Ky/M)G{Ay^N/M) ~ kA. 
Therefore, 

K[K'M-i^-^n)\nM-E^] 
■I - Ej} ^ ' 

3.3 Gyration radius 

The components of the gyration radius of a multi-block copolymer can be treated again by making 
use of the mapping on the homopolymer problem given by eqs H] and [3 In doing so the mapping 
looks as follows: 



u 



a — > aM 

N ^ ^f " 



2M 

Thus the gyration radius component in direction perpendicular to the surface becomes 

7e,x = aiV'^ (^A (^^y j (29) 

In the strong adsorption limit A^jN/M ^ 1 and 7^^ ~ aA^''^'^' M" , which yields 

2 30 

In a similar manner, the gyration radius component parallel to the surface has the form 

^.11 = ^^'^ ^.11 ( A (§) 1 (31) 



which in the hmit A-^N/M ^ 1 results in 

Like in the homopolymer case, one can define a blob length g^^ ~ {T^l/o)^^" ~ A^^/''^ M which 
in the strong adsorption limit, A > 1, approaches the block length, ^fefr — ^5 as it should be. 



Also in the limit of strong adsorption, A^jN/M ^ 1, the ratio 



/ AT \ 1^2 



- 



(33) 



\ M J \g,sj 
leads to the correct scaling in terms of number of blobs. 

4 Random copolymer adsorption 

Consider a random copolymer which is built up of Np A-type and i?-type monomers. The sam- 
pled AiJ-sequences are frozen (i.e. a distinct sample does not change during the measurement) which 
corresponds to quenched disorder. The binary variable a specifies the arrangement of monomers 
along the chain, so that a = 1, if the monomer is of A-type (A- monomers attract to the surface) 
and (7 = otherwise (i.e. in case of neutral 5-monomers). Let the fraction of attractive monomers 
(i.e., the composition) be p = Np/N and the fraction of neutral ones be 1 — p = N^/N . We assume 
that the statistics of sequences is governed by the Bernoulli distribution, i.e., the corresponding 
distribution function looks like: 

P {a] = p5{l - a) + {I - p)5{a) (34) 

This distribution is a special case of the more general Markovian copolymers^" when the " chemical 
correlation length" goes to zero. Two statistical moments which correspond to the distribution 
eq [31] are 

(a) = p 

{e') ^ {[a-{a)f)=p{l-p) (35) 



4.1 How does the critical Cc depend on the composition p? 

The adsorption of a random copolymer on a homogeneous surface has been studied by Whittington 
et al.^^'^^ within the framework of the annealed disorder approximation. Physically this means 
that during the measurements the chain touches the substrate at random in such a way that, as a 
matter of fact, one samples all possible distributions of monomers sequences along the backbone of 
the macromolecule. Following this assumption, let c'j^i'n) be the number of polymer configurations 
such that n units have contact with the surface simultaneously. The percentage of A-monomers 
(composition) is denoted by p. In the annealed approximation one then averages the partition 
function over the disorder distribution, i.e., 

N n 



n=l n„=0 ^ P 



N N 



5^4(n) K + 1 -P]" = E^^H (36) 



n=l n=l 



where e^g is the attraction energy of an effective homopolymer. From eq [36] one can see that the 
annealed problem is reduced to that of a homopolymer where the effective attractive energy is 
defined as 

e,\ = ln[pe^ + l-p] (37) 

We know that at the critical point the homopolymer attraction energy, e^^ = ej!, is model dependent. 
Then the critical attraction energy e = of a random copolymer reads 

e^c = In \'-^RA±P^] > (38) 

where the composition < p < I. Atp— >0e^— s>oo whereas at p = 1 = e^. The relationship in 
eq[38]has been recently found to be confirmed by Monte Carlo simulations.^^ 

5 Simulation Methods 



To check the theoretical predictions mentioned in the previous sections we have performed Monte 
Carlo simulations and investigated the adsorption of a homopolymer, multi-block copolymers, and 
random copolvmers on flat surfaces. Two coarse-erained models, the bead sprinE model and the 



simple cubic lattice model, Figure [21 are used, and two different Monte Carlo algorithms, the 
Metropolis algorithm (MA) and pruned-enriched Rosenbluth method (PERM), are applied to the 
two models, respectively. 




Figure 2: Schematic representation of a grafted chain close to criticality. (a) Snapshot of a chain 
with length = 128 from the MA model and block size M = 2; (b) A^ = 2048 with M = 8 from 
the PERM simulation. 



5.1 Off- lattice bead spring model with MA 

We have used a coarse grained off-lattice bead spring model^ to describe the polymer chains. Our 
system consists of a single chain tethered at one end to a flat structureless surface. There are two 
kinds of monomers: "A" and "B", of which only the "A" type feels an attraction to the surface. 
The surface interaction of the "A" type monomers is described by a square well potential U^iz) = e 
for z < 6 and Uw{z) = otherwise. Here e/ksT is varied from 0.6 to 3.6. The effective bonded 
interaction is described by the FENE (finitely extensible nonlinear elastic) potential. 

l-lo 



Ufene = -K{1 - lo) In 



^max ^0 



(39) 



with K = 20, Imax = 1, ^0 = 0.7, Imin = 0.4 



The nonbonded interactions are described by the Morse potential. 

— — = exp(-2a(r - rmin)) - 2exp(-a(r - r^™)) (40) 
Cm 

with a = 24, Tmin = 0.8, eM/ksT = 1. 

We use periodic boundary conditions in the x — y directions and impenetrable walls in the z 
direction. We have studied polvmer chains of lengths 32, 64, 128, 256 and 512. We have also 



studied homopolymer chains and random copolymers (with a fraction of attractive monomers, 
p = 0.25, 0.5, 0.75). The size of the box was 64 x 64 x 64 in aU cases except for the 512 chains 
where we used a larger box size of 128 x 128 x 128. The standard Metropolis algorithm was 
employed to govern the moves with self avoidance automatically incorporated in the potentials. In 
each Monte Carlo update, a monomer was chosen at random and a random displacement attempted 
with Ax, Ay, Az chosen uniformly from the interval —0.5 < Ax, Ay, Az < 0.5. The transition 
probability for the attempted move was calculated from the change AU oi the potential energies 
before and after the move as W = exp{—AU/kBT). As for standard Metropolis algorithm, the 
attempted move was accepted if W exceeds a random number uniformly distributed in the interval 



5.2 Coarse-grained lattice model with PERM 

The adsorption of AB block copolymer with one end (monomer A) grafted to a plane impenetrable 
surface and with only monomers A attractive to the surface are described by SAWs oi N — 1 steps on 
a simple cubic lattice with restriction z > 0. There is an attractive interaction between monomers 
A and the wall. The partition sum now is written as 



where A}sf{Ns) is the number of configurations of SAWs with N steps having Ng sites on the wall, 
and q = e^/*^*^ {ksT = 1 hereafter) is the Boltzmann factor, e > is the attractive energy between 
the monomer A and the wall. As g — > 1, there is no attraction between the monomer A and the wall. 
On the other hand it becomes clear that any copolymer will collapse onto the wall, if q becomes 
sufficiently large. Therefore we expect a phase transition from a grafted but otherwise detached to 
an adsorbed phase, similar to the transition observed also for homopolymers. 

For our simulations, we use the pruned-enriched Rosenbluth method (PERM)^^ which is a biased 
chain growth algorithm with resampling ("population control") and depth- first implementation. 
Polymer chains are built like random walks by adding one monomer at each step. Thus the total 
weight of a configuration for a polymer consisting of N monomers is a product of those weight 
gains at each step, i.e. Wn — H^o^Wj. As in any such algorithm, there is a wide range of possible 



[0,1]. 




(41) 



distributions of sampling, we have the freedom to give a bias at each step while the chain grows, and 
the bias is corrected by means of giving a weight to each sample configuration, namely, Wi Wi/pi 
where Pi is the probability for putting the monomer at step i. In order to suppress the fluctuations 
of weights as the chain is growing, the population control is done by "pruning" configurations with 
too low weight and "enriching" the sample with copies of high- weight configurations. Therefore, 
two thresholds are introduced here, = c+Z„ and W~ = c~Zn, where Zn = jf- J2config from 
the Mn trail configuration is the current estimate of partition sum at the n — 1 step, and c~ are 
constants of order unity and jc' ~ 10. In order to compare with the results obtained by the first 
method, we simulate homopolymers of length = 2048 and multi-block copolymers with block 
size M = 2^^ , = 0, 1, 2, ■ ■ ■ , 9. The number of monomers is increased to = 8192 as the block 
size increases. Also random copolymers of A^ = 2048 monomers with composition p = 0.125, 0.25, 
0.50, and 0.75 are sampled. 

6 Simulation Results 

6.1 Determination of the critical point of adsorption 

The determination of the critical adsorption point (CAP) is essential for testing the scaling results 
and for comparison with theory. In this work we determine the CAP from the analysis of several 
quantities: the order parameter /, the variance of the number of adsorbed monomers, C^, and the 
gyration radius Rg. These methods are described as follows: 

6.1.1 CAP from the order parameter 

From the plots of the order parameter / against the adsorption energy e for chains of different length 
A^ we determine the CAP as the point where the tangent taken at the inflection point of the order 
parameter curve intersects the horizontal axis e. Results are shown in Figure [3] for homopolymers 
and in Figure H] for multi-block copolymer with block size M = 2. In Figure [3^ and UK data is 
obtained by MA method in our off-lattice model, while in Figures |3b and HJd the data is obtained 
by PERM for self-avoiding chains on a cubic lattice. Evidently, in both cases the order parameter 
/ increases with growing strength of the substrate potential e. Thus the polymer chain undergoes 




Figure 3: The order parameter / against the adsorption energy e for homopolymers of different 
chain lengths N. The value of the CAP e'^lN) for — > cx3 is extrapolated from the log- log plot of 
e^(A^) versus as shown in the insert. In the thermodynamic limit (a) ~ 1.716 (MA off-lattice 
model), (b) ^ 0.284 (PERM on a cubic lattice). 

a transition from a grafted, but otherwise detached state, to an adsorbed state whereby the chain 
lies flat on the surface plane - see Figure [2b. The transition region narrows down as A^ increases, 
which is in good agreement with the scaling prediction of /, eq [161 iii all cases. In the inset of 
Figures [3] and Figure [U we see that the critical point t^{N) for homopolymers of chain length A^ 
as well as the critical points e^'\N) for multi-block copolymers of chain length A^ with M = 1, 
M = 2, M = 4, M = 8, and M = 16, gradually increase as A^ — oo. By extrapolating the 
data to 1/A^ = 0, one obtains the CAP values in the thermodynamic limit. Results for by MA 
and by PERM are listed in Tables [1] and [2l We should point out here that the simulation with 

] 



Table 1: MA 



M/N 


64 


128 


256 


512 


oo, / 


Rg 


1 


2.47(3) 


2.58(3) 


2.63(3) 


2.63(3) 


2.672(30) 


2.65(3) 


2 


2.32(3) 


2.44(3) 


2.47(3) 


2.48(3) 


2.52(2) 


2.52(3) 


4 


2.13(3) 


2.260(3) 


2.29(3) 


2.29(3) 


2.34(2) 


2.30(4) 


8 


1.93(3) 


2.08(3) 


2.12(3) 


2.14(3) 


2.19(3) 


2.06(4) 


16 


1.76(3) 


1.93(3) 


2.00(3) 


2.01(3) 


2.06(3) 


1.95(4) 


P/N 














1.0 


1.62(2) 


1.66(2) 


1.701(20) 


1.698(25) 


1.716(20) 


1.718(20) 


0.75 


1.83(2) 


1.89(2) 


1.92(2) 


1.946(20) 


1.95(3) 


1.95(3) 


0.50 


2.21(2) 


2.25(2) 


2.29(2) 


2.32(2) 


2.33(2) 


2.38(5) 


0.25 


2.81(4) 


2.97(4) 


2.98(4) 


3.02(4) 


3.05(5) 


2.91(6) 




Figure 4: The order parameter / plotted as a function of attractive energy e for copolymers with 
block size M = 2. The extrapolation plots for €c{N) versus 1/M for block sizes M = 1, 2, 4, 8, 16, 
and for the homopolymer, plotted versus are shown in the insert, (a) - the MA model, (b) - 
PERM. 



Table 2: PERM 



M/N 


64 


128 


256 


512 


1024 


oo, f 


Rg 


1 


0.337(9) 


0.457(5) 


0.505(5) 


0.535(3) 


0.548(2) 


0.560(2) 


0.568(6) 


2 


0.322(4) 


0.438(4) 


0.486(3) 


0.516(2) 


0.536(2) 


0.545(8) 


0.556(3) 
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0.296(7) 


0.411(4) 


0.465(3) 


0.489(3) 


0.511(2) 


0.520(4) 


0.524(2) 


8 


0.368(4) 


0.422(4) 
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0.464(3) 
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0.480(2) 


0.478(3) 
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0.411(2) 


0.426(3) 


0.432(2) 


0.441(2) 


0.437(4) 
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the MA model requires considerable computational effort for > 512, therefore, with the PERM 
method we confine ourselves to chain lengths not large than A^ = 2048 (Figure [3b), which are 
considerably shorter than feasible.^ Nevertheless, our estimate of the CAP e'l = 0.285(3) is in good 
agreement with previous results^ (within the error bars) although corrections to scaling have not 
been considered here. 



6.1.2 Prom the variance of the order parameter: 

In a computer simulation one usually computes the variance of the order parameter, A/, which 
yields some important thermodynamic quantities like isothermal compressibility, and/or specific 
heat, via the fiuctuation relations. 

N'Af = (iVf) - (iV,)2, (42) 

At the CAP A/ has a maximum which becomes larger and narrower as one approaches the thermo- 
dynamic limit, N ^ oo. In Figure [5^ this is shown for the PERM model along with an extrapolation 
of the CAP e^(A^) for chains of length - see inset - which for N oo becomes a straight line 
in agreement with eq [161 It becomes also evident from Figure [5]d that the alternative method of 
using the position of the maximum of the specific heat, Cy = [kBT'^)"^ {{U"^) ~ (U)"^) from the 
fiuctuations of the internal energy, U = eNg, does not give satisfactory results due to the rather fiat 
shape of the maximum. This behavior is not surprising, if one recalls that the critical exponent a 
describing the divergence of Cy at the CAP, i.e., for /t — > 0, according to Cy oc see eq[T5l is 
given by a = 2 — (p^^ ^ 0.^^ It has been show earlier, however, that one can still use specific heat 
data to determine the CAP if, instead of the position of the maximum, one examines the common 
intersection point of Cy vs e. In our simulation this yields again = 0.284 - cf. Table O 

6.1.3 Prom the components of Rg 

According to eqs [6], [9l and [131 one should expect that all curves of Rl±/Rgp for different chain 
length intersect at a fixed point which gives the CAP in the limit of A^ — > oo. In Figure [6l we 
illustrate this method by plotting the ratio R'^j_/R'^^^ vs e for copolymers with block size M = 2. For 
both methods, MA and PERM, the curves for different A^ intersect nearly at a single intersection 
point, however, as before, the CAP determined by MA (see Figure [6^) is less accurate than the 
results given by PERM (see Figure [6b). The CAPs obtained from this method, e^^^^ = 2.52(3) by 
MA and ef^^^ = 0.556(4) by PERM are consistent with the estimates from the order parameter 
method where ef =2 = 2.521(20) by MA and ef =2 = 0.546(8) by PERM. The CAPs e^M) for 
homopolymers, multi-block copolymers with different block size M, and for random copolymers are 
listed in Table [Tl and [2l 
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Figure 5: Results for (a) the variance of the order parameter / multiphed by in the case of a 
homopolymer. The inset shows an extrapolation of the CAP e^(A^) for — * which converges 
to the value for an infinite chain, ej? = 0.284 - cf. Table [21 (b) the specific heat per monomer, 
Cy/N, which plotted as a function of e for homopolymers of different chain length N."^"^ 



6.2 Scaling behavior 

From the data for the CAP one may check the value of the crossover exponent (j) = 0.50 by plotting 
the order parameter / vs. A^ (eq[3]). This is illustrated in Figure [7] as a double logarithmic plot of / 
vs. A^ for the case of M = 1, i.e., regular alternating polymers. Figure [7] demonstrates clearly that 
the slope of the f vs N curves in logarithmic coordinates is equal to 1 — = —0.5 only in those 
cases where the strength of the substrate potential equals the CAP value ec, in agreement with the 
relation / oc N'f'"^. As in the case of homopolymers (eq[3]). Figure [7^ shows that in the strongly 
adsorbed regime (e = 3.40) above the CAP the order parameter / oc A^*^ {all monomers stick). In 
contrast, far below the CAP, only the anchoring monomer is attached to the substrate, / ~ A^~^, as 
in the asymptotic limit A^ oo of homopolymers. This is observed for e = 0.60 for the alternating 
chains (M = 1). In Figure [7|3, where the statistical precision and the chain lengths involved are 
much higher, one may see that for large A^ the curves which are slightly above, e = 0.571, or below, 
e = 0.560, the CAP at e = 0.568 - cf. Figure 0- display slopes which differ slightly from —0.5 and 
thus considerably narrow the interval of critical behavor. 

In Figure [8] we present the results for the components of the mean square gyration radius, R^h 



and R^±, in scaled form in terms of the parameter kN'^ for regular block-copolymers with block size 
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Figure 6: The ratio of -R^xZ-R^y plotted as a function of e for copolymers with block size M = 2. 
The critical point is determined by the intersection of all curves which are found to be at (a) 
ei*'^=^^ ^ 2.52 (MA) and (b) ef^=^^ ^ 0.556 (PERM). 



M = 1 and M = 8. Generally, one observes a good agreement with the predictions of Section [21 
especially concerning the data obtained by PERM - Figure [Hb, d. Considerable deviations from the 
expected scaling behavior are observed only in Figure [8]d where the effective segment of a diblock 
with M = 8 is comparatively large for the simulated chain lengths iV < 512, meaning effective chain 
lengths of N^ff = N/16 < 32 which are definitely too short for a well pronounced scaling behavior 
to be demonstrated. With the much longer chains, < 2048, sampled by PERM and shown in 
Figure [Hll, this problem is absent. 

6.3 Phase diagram of multi-block copolymer adsorption 

Using the values for the CAP, given in Tabled], one may construct a phase diagram showing the 
relative increase of the critical potential ec(M) compared to that of a homopolymer against (inverse) 
block size M. This is one of the central results of the present study. In Figure [9] one may see that 
the line of critical points, defining the region of adsorption, for both models is a steadily growing 
function of the inverse block size M^^. Evidently, the theoretical result, eq[25], appears to be in 
good qualitative agreement with simulation data for the different models. As far as eq [25] comes as 
a result of scaling analysis, it can be verified only up to a factor of proportionality. As mentioned 
in Section [2?n the CAP of a homopolymer, e^, is of the same order as that of the "renormalized" 
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Figure 7: Log-log plot of the order parameter / vs N for block copolymers with block size (M = 1). 
The value of e for each curve is given in the legend while the slope is also indicated. One may 
readily check that the straight lines with slope 0.5 correspond to the respective values of ec in both 
models, (a) MA and (b) PERM. 



chain consisting of diblocks, E^. Thus from a fit of the data points with the expression eg 1251 one 
may actually determine E^. So in the MA model one gets E'^ = 3.306 and for PERM E'^ = 1.254, 
that is, one gets values which are two to four times larger than the respective CAP values of a 
homopolymer in both models. 

6.4 Random Copolymers 

In this section we examine the adsorption transition of random copolymers with quenched disorder 
and average percentage p of the A monomers. In addition to testing the scaling behavior, we also 
check to what extent one may employ the theory developed within approximation of "annealed 
disorder" for the description of the CAP properties. We performed Monte Carlo simulations for 
heterogeneous random copolymers of chains lengths 32, 64, 128, 256 and 512 (MA) and for 64 < 
< 2048 (PERM) with different fraction of attractive monomers {p = 0.125,0.25,0.50 and 0.75). 
It has been pointed out earlier^^'^^ that the crossover exponent stays the same, = 0.5, also in 
the case of random copolymers. Both simulation methods used in the present study demonstrate 
this in Figure [10] where qualitatively the observed picture is similar to that of Figure [7] - small 
deviations in the attraction potential e, which was used when sampling the values of the order 
parameter /, manifest themselves in significant changes of the log-log slope 1 — from the expected 
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Figure 8: Log-log plots of Rln/N'^" and R^JN^^ vs nN''' with z/ = 0.588 and 1/3 = 3/4. The 
straight lines indicate the asymptotic behaviour of the scaling functions given by eq [6] and [9j (a) 
and (b) represent results for regular multi-block copolymers with block sizes M = 1 and M = 8, 
respectively, and are obtained by MA; (c) and (d) - similar results but obtained by PERM. 



value of —0.5. In Figure [TT] we demonstrate that the scaling of the mena square gyration radius 
components, which we discussed before with regard to the multiblock copolymers, holds also for 
random copolymers with different composition p. Again the value of = 0.5 gives best scaling 
results. Thus it turns out that the composition affects only the value of the CAP e^. 

In Figure fT2] we present a plot of the critical point of adsorption against the fraction of attractive 
monomers. The full line corresponds to the theoretical prediction, eq[38l Given that there are no 
fitting parameters in this equation, one finds a very good agreement between theoretical predictions 
and simulation results as well as with very recent simulation results^^ which demonstrates the 




Figure 9: = (e*^ — e^)/t^ plotted vs 1/M for multi-block copolymers with various values of M. 
The critical points of adsorption for homopolymers are (a) = 1.716 (MA) and (b) = 0.285 

1 /2 

(PERM). The curves give the best fit of eq ESI k oc ft^MMlt^^ 



the adsorption of random copolymers can be properly described within the scope of the annealed 
approximation. Figure [T2] also indicates that this approximation breaks down for chains which 
are not random^^ - at 50% composition the CAPs of regular block copolymers are clearly off the 
theoretical prediction, eq[38l As far as polymer adsorption is greatly facilitated by the formation of 
trains of monomers on the substrate/^ the larger the block size M, the lower the respective CAP 
under the line, eq[38l No monomer trains are possible in the case of alternating chains which 
results in an e^^^^ > e^. Thus from the position of the CAPs on Figure [T2] one may conclude that 
the mean length of an A-train on the substrate at p = 0.5 is close to four. 

7 Concluding remarks 

The main focus of the present investigation has been aimed at the adsorption transition of random 
and regular multiblock copolymers on a rigid substrate. We have used two different models to 
establish an unambiguous picture of the adsorption transition and to test scaling predictions at 
criticality. The first one is an off-lattice coarse-grained bead-spring model of polymer chains which 
interact with a structureless surface by means of a contact potential, once an A-monomer comes 
close enough to be captured by the adsorption potential. The second one deals with SAW on a 
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Figure 10: The same as in Figure [7] but for random copolymers with the composition p = 0.5 (a) 
MA where = 2.33 and (b) PERM with eP = 0.507. 

cubic lattice by the pruned-enriched Rosenbluth method (PERM) which is very efficient, especially 
for very long polymer chains, and provides high accuracy of the simulation results at criticality. 
Notwithstanding their basic difference, both methods suggest a consistent picture of the adsorption 
of copolymers on a rigid substrate and conffim the theoretical predictions even though the particular 
numeric values of the critical adsorption potential (CAP) are model-specific and differ considerably. 

As a central result of the present work, one should point out the phase diagram of regular 
multiblock adsorption which gives the increase of the critical adsorption potential e*^ with decreasing 
length M of the adsorbing blocks. For very large block length, M^^ — > 0, we find that the CAP 
approaches systematically that of a homogeneous polymer. We demonstrate also that the phase 
diagram, derived from computer experiment within the framework of two different models, agrees 
well with the theoretical prediction based on scaling considerations. 

The phase diagram for random copolymers with quenched disorder which gives the change in 
the critical adsorption potential, e^, with changing percentage of the sticking A-monomers, p, is 
also determined from extensive computer simulations carried out with the two models. We observe 
perfect agreement with the theoretically predicted result which has been derived by treating the 
adsorption transition in terms of the "annealed disorder" approximation. 

We show that a consistent picture of how some basic polymer chain properties of interest such 
as the gyration radius components perpendicular and parallel to the substrate, or the fraction of 





Figure 11: The same as in Figure [8] but for random copolymers at different composition p. 

adsorbed monomers at criticality, scale when a chain undergoes an adsorption transition appears 
regardless of the particular simulation approach. An important conclusion thereby concerns the 
value of the universal crossover exponent (p = 0.5 which is found to remain unchanged, regardless 
whether homo-, regular multiblock-, or random polymers are concerned. 
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Figure 12: The CAP, e^, plotted vs the composition p for random copolymers. The curves give 
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